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ABSTRACT 

For any assumed standard stellar Initial Mass Function, the Sloan Digital Sky Survey 
(SDSS) gives a precise determination of the abundance of galaxies as a function of 
their stellar mass over the full stellar mass range 1O 8 M < M* < 10 12 M Q . Within the 
concordance ACDM cosmology, the Millennium simulations give precise halo abun- 
dances as a function of mass and redshift for all halos within which galaxies can form. 
Under the plausible hypothesis that the stellar mass of a galaxy is an increasing func- 
tion of the maximum mass ever attained by its halo, these results combine to give 
halo mass as a function of stellar mass. The result agrees quite well with observational 
estimates of mean halo mass as a function of stellar mass from stacking analyses 
of the gravitational lensing signal and the satellite dynamics of SDSS galaxies. For 
M* ~ 5.5 x 10 10 A/ Q , the stellar mass usually assumed for the Milky Way, the implied 
halo mass is ~ 2 x 1O 12 M , consistent with most recent direct estimates and inferences 
from the MW/M31 Timing Argument. The fraction of the baryons associated with 
each halo which are present as stars in its central galaxy reaches a maximum of 20% 
at masses somewhat below that of the Milky Way, and falls rapidly at both higher 
and lower masses. These conversion efficiencies are lower than in almost all recent 
high-resolution simulations of galaxy formation, showing that these are not yet viable 
models for the formation of typical members of the galaxy population. When inserted 
in the Millennium-II Simulation, our derived relation between stellar mass and halo 
mass predicts a stellar mass autocorrelation function in excellent agreement with that 
measured directly in the SDSS. The implied Tully-Fisher relation also appears consis- 
tent with observation, suggesting that galaxy luminosity functions and Tully-Fisher 
relations can be reproduced simultaneously in a ACDM cosmology. 

Key words: cosmology: theory - cosmology: dark matter mass function - galaxies: 
luminosity function, stellar mass function - galaxies: haloes - cosmology: large-scale 
structure of Universe 



1 INTRODUCTION 

It has long been known that baryons do not dominate the 
mass in the Universe, and recent observations of microwave 
background fluctuations have demonstrated that only about 
15% of cosmic matter is in the form of baryons. The re- 
maining 85% is apparently made up of some as yet unidenti- 
fied, weakly i nteracting, non-baryo nic particle, so-called cold 
dark matter |Komatsu et al.|[2009h . Together these two mat- 
ter components account for only a quarter of the current 
energy density of the Universe; the rest comes from a near- 
uniform dark energy field whi ch drives the obser ved acceler- 
ation of the cosmic expansion |Astier et alj|200d h Nonlinear 
structure formation in this concordance ACDM cosmology 
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proceeds through gravitationally driven hierarchical collapse 
and aggregation. Galaxies form by the cooling and conden- 
sation of gas at the centres of an evolving population of dark 
matter halos, a s orig inally set out in a different context by 
IWhite fc Reesl l|l97Sft . 

The abundance of dark matter halos as a function 
of their mass can be pr edicted by the simple theory of 
iPress fc Schechterl l|l974 ) and its extensions, and is ex- 
tremely broad. For example, halos with mass exceeding 
1O 14 M0 are predicted to contain roughly 10% of all dark 
matter in today's Universe; halos with mass below 1O" 5 M 
should contain another 10%; and the remaining 80% is 
distributed over the intervening 19 orders o f magnitude, 
with the median point ne ar 10 10 M Q (e.g. iMo fc Whitel 
l2002l : lAngulo fc Whitd[2009l h This behaviour has been con- 
firmed by detailed numerical simulations of cosmic evolution 
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over the limited mass range accessible to them (currently 
> 10 9 Mq; see below). The observed distribution of stars 
among galaxies is much more confined, however, with 10% 
of stars in galaxies with stellar mass above 1.3 x 10 n M o , 
10% in galaxies with stellar mass below 5 x 10 9 Mp, and a me- 
dian point near 4x 1O 1O M (|Li fc White!l2009l '). This implies 
that the baryons are converted into stars with very differ- 
ent efficiencies in halos of different mass. Star formation is 
most efficient at the centres of halos of typical galaxies such 
as the Milky Way, and is substantially less efficient at the 
centres of much more mass i ve or much le ss massive halos 
|Navarro fc Steinmet3l200Ch . iLarsonl 11974 ) noted that star 
formation would likely be inefficient in objects with low es- 
cape velocitie s because of the ease with which winds can 
expel gas, and I White fc R ces (1978) invoked this process to 
explain the relatively small fraction of all stars which end up 
in low-mass galaxies. The inefficiency of star formation at 
the centre of massive halo s is related to the we ll known "cool- 
ing flow paradox" (e.g. Fabian et al.l [2001). Recent work 
suggests tha t it may result from feedback from the central 
radio AGN jTabor fc Binnevlll993l; ICiotti fc Ostrikerll200ll ; 
iBirzan et al.[|2004l ; ICroton et al.ll2006l ; iBower et alJ^OOrj ) 

There are several ways to study how galaxies popu- 
late dark matter halos. The most straightforward is di- 
rect simulation of galaxy formation in its cosmological con- 
text. An A-body treatment of the evolution of the dark 
matter comp onent is combined w it h either a hydrody 
namical (e.g. ICen fc Ostrikerl |2000|; ISpringel fc Hernquist 



2003 1; iKeres et al]|2005l ; iPfrommer et al.ll2006l; ISiiacki et al 
20071 ') or a semi-analytic (e.g. Kauffma nn et al 
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Springel et ai]|200ll ; lHatton et aT]|2003l ; ISpringel et ai1l2005l ; 
Kang et al.ll2005r ) treatment of baryonic evolution. The ad- 



vantage of these methods is that they track galaxy and dark 
matter halo evolution across cosmic time in a physically con- 
sistent way, providing positions, velocities, star formation 
histories and other physical properties for the galaxy popu- 
lations of interest. In recent years both techniques have had 
considerable success in reproducing observations. Hydrody- 
namical simulations provide a much better description of 
diffuse gas processes, but are relatively inflexible, typically 
producing galaxy populations which are at best a rough 
match to observation. The greater flexibility and speed of 
semi-analytic methods allows much better reproduction of 
observational data at the expense of a schematic treatment 
of diffuse gas physics. Both schemes implement simplified 
and highly uncertain recipes to treat star and black hole 
formation and related feedback processes. For both, these 
entrain considerable uncertainty about whether the physics 
of galaxy formation is reliably represented. 

Even better fits to the observed luminosity, colour and 
clustering distributions of galaxies can be obtained by giv- 
ing up on any attempt to represent formation physics, and 
instead using simple statistical models with adjustable pa- 
rameters to populate dark matter halos with galaxies. By 
adjusting model assumptions and their associated param- 
eters, this halo occupation distribution (HOD) approach 
is able to match observed statistics like galaxy luminosity 
and correlation functions as a function of luminosity and 
colour with remarkable accuracy within the concordance 
ACDM cosmology (e.g . ICooravl I2OO2I ; iBerlind fc Weinberg! 
l2002i ; lYang et al.l 2003h . Their disadvantage is that they use 
no information about the evolution of a system when pop- 



ulating it with galaxies, and this formation history may 
have sig nificant influence on the properties of the galax- 
ies (e.g. iGao et alj [20051 ; ICroton et all 120071 ). Several au- 
thors have explored schemes intermediate between direct 
simulation and HOD modelling in an attempt to retain 
some of the advantages of each (e.g. IWang et al.ll200r3 . 120071 ; 
IConrov fc Wechsled[2009h . 



An alternative is to link galaxies to their halos/subhalos 
by matching observed galaxy luminosity functions to sim- 
ulated halo mass functions assuming a unique and mono- 
tonic relation between galaxy lu minosity and hal o mass . 
This method was proposed by IVale fc Ostrikerl {2004) 
and then extended by several groups to consider a va- 
riety of propert ies both for the g a laxies and for their 
halos/subhalos JConrov et all 120061; IShankar et al.l 120061 ; 
IConrov et"aT]|2007l; iBaldrv et al.ll2008l ; iMoster et alj|2009[ ). 
In particular. IMoster et al.l |2009l ) show how the method can 
be extended to allow a scatter in the properties of galaxies 
associated with halos or subhalos of given mass (see also 
I Wetzel fc Whitell2009l ). 

Here we follow this latter approach. We combine a pre- 
cise stellar mass function based on the full s pectroscopic 
dataset of the most recent SDSS data release (|Li fc White! 
2009) with a precise halo/subhalo m ass function obtained 
from the Millennium Simulation CMS; ISpringel et al.ll2005h 
and the higher resolution M illennium-II Simulation (MS-II; 
iBovlan-Kolchin et al.l [20091 '). This yields a much more ac- 
curate relation between galaxy stellar mass and dark mat- 
ter halo mass than could be derived from earlier data. We 
compare this relation with direct observational estimates of 
the mean mass of halos surrounding galaxies of given stellar 
mass inferred from gravitational lensing and satellite galaxy 
dynamics data. We also compare the halo masses predicted 
for the Milky Way and other Local Group galaxies with esti- 
mates derived from dynamical data. As a further consistency 
test we populate halos/subhalos in the MS and MS-II with 
galaxies of stellar mass chosen according to our relation, and 
we compare the stellar mass correl ation function of th e re- 
sult with the SDSS measurement of iLi fc White] (|2009l V We 
derive the star formation efficiencies implied as a function of 
halo mass by our relation, and we compare them to the effi- 
ciencies in published hydrodynamical simulations of galaxy 
formation in the ACDM cosmology. Finally, we revisit the 
issue of whether the abundances and circular velocities of 
galax ies can be fit s imultaneously in a ACDM cosmology 
fe.g. lCole eTal]|2000! ). 

In the next section we briefly describe the two Millen- 
nium simulations, and the definitions of halo mass that we 
adopt for the rest of this paper. The relation between galaxy 
stellar mass and dark matter halo mass is derived at the be- 
ginning of our results section, Sec. [3] and in later subsections 
we compare it with direct observational determinations, we 
give predictions for the halo masses of Local Group galax- 
ies, we compare the implied stellar mass correlation function 
with that measured in SDSS, and we revisit the problem of 
simultaneously reproducing the luminosity function and the 
Tully-Fisher relation in a hierarchical cosmology. Conclu- 
sions and a discussion of our results are presented in Sec. [4] 
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2 DARK MATTER HALOS 

Both the Millennium Simulation and the Millennium-II 
Simulation adopt the concordance ACDM cosmology with 
parameters c hosen to agree wit h a combined analysis of 
the 2 dFGRS dColless et aUkoOll ) and the first-year WMAP 
data (|Spergel et al. 20031 ). The parameters are f2 m = 0.25, 
Q b = 0.045, H = 73km S - 1 Mpc- 1 , ft A = 0.75, n = 1, 
and <rg = 0.9. These parameters are only marginally con- 
sistent with an alysis of the latest W MAP and large-scale 
structure data (|Komatsu et al.|[2009l ) but the differences are 
too small to significantly affect the analysis of this paper. 
Both simulations use 2160 3 dark matter particles to trace 
evolution from z ~ 127 to z ~ 0. They were carried out in 
cubic volumes with periodic boundary conditions and sides 
of length 685 Mpc and 137 Mpc for the MS-I and MS-II, re- 
spectively, corresponding to particle masses of 1.2 x 10 9 M o 
and 9.5 x 10 6 Mq. The large volume of the MS enables one to 
study even the cD galaxies of rare and massive clusters with 
good statistical power, while the excellent mass resolution 
of the MS-II can resolve the dark matter halos predicted to 
host even the faintest known dwarf galaxies. 

For each of the output dumps, friends-of-friends (FOF) 
groups are identified in each simulation by linking together 
particles separ ated by less than 0.2 of the mean interparti- 
cle separation jDavis et al.lll985h . The SUBFIND algorithm 



( Spr ingel et al 



200 ll ) was then applied to each FOF group 



in order to split it into a set of disjoint, self-bound subha- 
los, which represent locally overdense and dynamically sta- 
ble subunits within the larger system. The main subhalo is 
defined as the most massive self-bound subunit of a FOF 
group and normally contains most of its mass. Merger trees 
were then built which link each subhalo present in a given 
dump to a unique descendent in the following dump. These 
allow us to track the formation history of every halo/subhalo 
pre sent at z = 0. We refe r read ers to ISpringel et alj (|2005T ) 
and lBovlan-Kolchin et all (|2009h for a more detailed descrip- 
tions of these simulations and post-processing procedures. 

In this work, we assume that both main subhalos and 
satellite subhalos have galaxies at their centres, and that 
the stellar masses of these galaxies are directly related to 
the maximum dark matter mass ever attained by the sub- 
halo during its evolution. We denote this mass by Mhaio- 
In practice this mass is usually the mass at z = for main 
subhalos, and the mass just prior to accretion for satellite 
subhalos. Semi-analytic simulations show that for satellite 
systems this latter mass is much more closely related to 
the stellar mass of the central galaxy than is the z — 
mass of the subhalo, because the latter has often been very 

3 



2001 



substantially reduced by tidal st ri pping dGao et al. 
Wang et alJ l2 006r,lFont et alJl2008h.lvale fc Ostrikerl l|2004h . 
Conrov et al.l (|2006t ) and Berrier et al.l (|2006t ) present gen- 



eral plausibility arguments for such an assumption rather 
than studying any specific galaxy formation model. We thus 
need to estimate the abundance of (sub)halos in the Millen- 
nium simulations as a function of this Mhaio ■ 

For each FOF group, we define the centre as the min- 
imum of the gravitational potential well and we define the 
virial radius, Rvir, as the radius that encloses a mean over- 
density of 200 times the critical value. The mass within R v i r 
is then defined as the virial mass: 
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Figure 1. Dark matter halo mass functions at z = where "halo" 
is defined to include both main subhalos and satellite subhalos. 
Halo mass. Mhaio, ls defined as the current virial mass for main 
subhalos and as the virial mass immediately prior to accretion for 
satellite subhalos. In both cases this is normally the maximum 
mass attained over the subhalo's history. Black triangles are for 
the MS and red dots are for the MS-II. Poisson error bars based 
on halo counts are shown for both simulations. 



M ha 



100 a 3 
— H (z)R vl 



(1) 



We define Mhaio for a main subhalo to be its current virial 
mass, and for a satellite subhalo to be its virial mass im- 
mediately prior to accretion onto a larger system, i.e. its 
virial mass immediately before it last switches from being a 
main subhalo to a satellite subhalo. Hereafter, we refer to 
both main subhalos and satellite subhalos as "halos" and we 
refer to Mhaio defined in this way as the "halo mass" . 

FigH] shows halo mass functions for the two Millennium 
simulations at z = 0. Black triangles refer to the MS and 
red dots to the MS-II. The two simulations agree well above 
1O 12 ' 3 M0 but below this threshold, the MS lies progressively 
below the MS-II. This is due to resolution effects, which set 
in at substantially h igher masses than for the FO F halo mass 
function in Fig. 9 of iBovlan-Kolchin et all (2009) because of 
the inclusion of satellite subhalos. These can fall below the 
resolution of the MS at z — yet still be relatively massive 
at the time of infall. In the following, we combine the part 
of the MS mass function with M ha io > 1.9 x 10 12 M Q with 
the part from the MS-II with M ha io < 1.9 x 1O 12 M in 
order to represent the overall dark halo mass function as 
well as possible. Based on the deviations between the two 
simulations visible in Fig[TJ we estimate that the resulting 
function should be accurate to better than about 10% from 
1O 1O M up to 10 15 A/ Q . This will turn out to cover the full 
halo mass range of relevance for real galaxies. 
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Figure 2. The stellar mass - dark matter halo mass relation. 
The solid curve is obtained by matching galaxy abundances from 
SDSS/DR7 to dark matter halo abundances from the combina- 
tion of the MS and the MS-II (FigQl. The dashed curve shows 
an extrapolation of this relation to stellar masses of 10 6 Mq and 
1O 12 M0 at the low- and high-mass ends, respectively. The bottom 
panel shows the ratio of stellar mass to halo mass as a function 
of halo mass. 



band "model" luminosities. The latter are generally thought 
to give a better estimate of the total luminosity of galax- 
ies. This conversion is discussed in detail in Appendix A, 
which also gives a modified version of the fitting formula of 
iLi fc Whitefl|200Sl ) which represents this "total stellar mass" 
function. The correction increases stellar masses by about 
9% on average. If we leave aside uncertainties in the stellar 
Initial Mass F unction then results in Appendix A and in the 
appendices of ILi fc White! (|2009l ) suggest that the remain- 
ing systematic uncertainty in the stellar mass functions are 
of order 10% in stellar mass. Purely statistical errors are 
much smaller than this. Since the abundances matched in 
Eq. [2] range over almost six orders of magnitude, such un- 
certainties have only a very small effect on the Mhaio - M, 
relation. 

Our relation between galaxy stellar mass and the dark 
matter halo mass is shown in the upper panel of Fig. [2] 
The solid curve uses SDSS/DR7 data over the stellar mass 
range from 10 8 ' 3 to 10 11,8 Mq, which corresponds to dark 
matter halo masses between 1O 1O8 M and 1O 14 ' 9 M . We 
extrapolate this relation to 10 6 Mq at the low-mass end and 
to 10 12 Mq at the high-mass end assuming constant slope, 
as indicated by the dashed extensions. Galaxies with mass 
around 10 6 Mq are expected to reside in dark matter halos 
with mass ~ lO lo Af0, where we expect errors in our abun- 
dance estimates still to be below 10%. At the high-mass end, 
the stellar mass of the central galaxy becomes very insensi- 
tive to its dark matter halo mass, indicating a suppression 
of star formation in the cores of halos more massive than 
~ 1O 13 M . 

If we a dopt the funct i onal f orm suggested bv lYang et al.l 

(|2003h and lMoster et all lj2009h . our derived relation can be 
approximated to high accuracy by 



3 GALAXY FORMATION EFFICIENCY 

3.1 Connecting Galaxies to Dark Matter Halos 

We connect dark halo mass Mhaio to the stellar mass of the 
central galaxy by assuming a one-to-one and monotonic rela- 
tionship between the two. In practice, if the number density 
of dark matter halos with mass exceeding Mhaio matches 
the number density of galaxies with stellar mass exceeding 



M. 

Mhaio 



n(> Mhaio) = n(> M»), 



(2) 



then we assume galaxies of stellar mass M, to reside at the 
centre of dark matter (sub)halos of mass Mhaio- 

To derive the relation between Mhaio and M* we need to 
combine the halo mass function of Fig. [TJ with an equally pre- 
cise observed stellar mass function for galaxies . We ta ke the 
recent measurement presented bv lLi fc White! ((2009) . This 
is based on a complete and uniform sample of almost half 
a million galaxies from the Sloan Digital Sky Survey data 
relea se 7 (SDSS/DR7) ijYork et all l2000l : lAbazaiian et all 
2008). This extends over almost four orders of magnitude 
in stellar mass (10 8 Mq — 10 11,7 Mq) with very small sta- 
tistical errors. The main residual uncertainty comes from 
possible systematic errors in the determination of stellar 
masses from the SDSS photometry. Here we convert from 
masses based on SDSS r-b and Petrosian luminosities, as 
used by ILi fc White] (|2009h . to masses based on SDSS r- 



( M h alo \- a ( Mhalo \ 

I Mo j V M J 



(3) 



where c = 0.129, M = 10 114 M Q , a = 0.926, /3 = 0.261 
and 7 = 2.440. Note that this formula has been fitted to 
our results over the halo mass range 10 10 ' 8 to 10 14 ' 9 Mq, 
corresponding to the solid curve in the upper panel of Fig. [21 

In the bottom panel of Fig. [2] we show the ratio of 
stellar mass to dark halo mass as a function of halo mass. 
This reaches a maximum in halos with Mhaio ~ 10 11 ' 8 Mq, 
slightly less massive than the halos which host L* galaxies. 
The peak value is around 3.5%. The ratio drops very rapidly 
towards both lower and higher halo masses: M,/ Mhaio < 
0.27% in dark matter halos with mass ~ 10 10,7 Mq and 
M» /Mhaio ~ 0.09% in clusters with ~ 10 14 8 M Q . (Note 
that in the latter case the stellar mass refers only to the 
central galaxy.) 

Semi-analytic models like that of De Lucia & Blaizot 
(2007; DLB07) produce curves very similar to those of Fig. 2 
but noticeably offset. This offset comes from several sources. 
As may be seen in Fig.l, the MS does not produce the cor- 
rect (sub)halo abundance below Mhaio = 10 12 A/q, so that a 
SAM based on the MS alone (like that of DLB07)is skewed 
as a result. At the moment, there are no semi-analytic mod- 
els tuned to work simultaneously on the MS and MS-II sim- 
ulations, though we intend to produce such models in the 
future (Guo et al., in preparation) . In addition , the D LB07 
models do not accurately fit the ILi fc White! (|2009T ) mass 
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Figure 3. Dark matter halo mass as a function of stellar mass. 
The thick black curve is the prediction from abundance matching 
assuming no dispersion in the relation between the two masses. 
Circles with error bars are weak lensing estimates of the mean 
halo mass of central gala xies as a function of their stellar mass 
jMandelbaum et al.ll2006f ). The error bars show the 95% confi- 
dence ranges. Triangles with lcr error bars show mean halo masses 
as a function of central galaxy stellar mass derived fro m the 
stacked kinematics of satellite galaxies jMore et al.ll2009l) . Red 
and green dashed curves show abundance matching predictions 
for mean halo mass as a function of galaxy stellar mass assuming 
dispersions of 0.1 and 0.2, respectively, in logM* at given halo 
mass. The dashed black curve is the satellite fraction as a func- 
tion of stellar mass, as labelled on the axis at the right-hand side 
of the plot. 



function (see their Fig. Bl) so this also introduces an appre- 
ciable offset in the M* - Mhaio relation. The scatter of this 
model around its own M* - Mhaio relation is, however, quite 
small and is comparable with the values we test below. 



3.2 Comparison to observed Mhaio - M 4 relation 

We now focus on galaxies with stellar mass between 10 10 Mq 
and 3 x 1O 11 M0 and show halo mass vs. stellar mass in 
Fig. [3] The solid curve is the prediction of our abundance 
matching as shown already in Fig. [2] The circles with er- 
ror bars show mean halo mass as a function of central 
galaxy stellar ma ss as obtained from gra vitational lens- 
ing measurements (jMandelbaum et al.ll2006l ). The data con- 
sisted of 351,507 galaxies from SDSS, including both early- 
and late- type galaxies in the mass range [0.7 x 10 10 Mq, 
4 x 10 Mq]. The error bars indicate 95% confidence inter- 
vals. The stellar masses used here are based on the pho- 
tometric prop erties of SDSS galaxies (see Appendix A), 
while those in iMandelbaum et all (|2006l ) were estimated by 
Kauffmann et al. I (|2003h based on the SDSS spectroscopy. 
We have shifted the latter by about 0.1 dex to match the 
form er. This conv e rsion factor was obtained as in Appendix 
A of iLi fc Whitel ([2009) using a large sample of galaxies 
for which both stellar mass estima tes are available (see also 
Fig. 17 of lBlanton fe Roweisl (2007)). The lensing and abun- 



dance matching results agree well for stellar masses between 
10 10 M o and 1O 1O - 8 M . At higher masses, the abundance 
matching halo masses are somewhat larger than those esti- 
mated from lensing. 

Triangles with error bars in Fig. [3] show mean halo mass 
as a function of central stellar mass as derived from the 
stacked kinemati cs of the satellites of a large sample of SDSS 
central galaxies |More et al.l 12009). Error bars here denote 
lcr uncertainties in the mean. To make a dire ct compari- 
son w ith our predictions, we have converted the lMore et al.l 
(2009) relation between halo mass and galaxy luminosity 
into a relation between halo mass and central stellar mass 
by using the SDSS data to estimate mean "model" stellar 
mass-to-Petrosian light ratio as a function of stellar mass. 
We have also accounted for our differing definitions of "halo 
mass" using an NFW profile of appropriate concentration. 
The resulting Mhaio - M, relation is consistent with that 
found from the lensing measurements, although its slope is 
slightly flatter. Again, the Mhaio inferred at high mass is 
somewhat lower than predicted by our abundance match- 
ing. 

We have studied whether this discrepancy could be due 
to a significant dispersion in the central stellar mass of ha- 
los of given dark matter mass. Assuming the dispersion in 
logM* to be gaussian and independent of Mhaio, rms values 
exceeding about 0.2 are excluded because they are incon- 
sistent wit h the steep high-m ass fall-off of the stellar mass 
function of ILi fc Whitel (|2009l ). 

In Fig. [3] we plot the mean halo masses predicted for 
distributions which match both the SDSS/DR7 stellar mass 
function and the halo mass function of Fig. [1] assuming dis- 
persions in logM* at given Mhaio of 0.1 (red curve) and 0.2 
(green curve); our standard relation (the thick black curve) 
corresponds to a dispersion of zero. Dispersions in the al- 
lowed range lead to a flatter slope of the Mhaio - M, rela- 
tion and to the prediction of a lower mean halo mass at high 
central stellar mass. The effects are quite weak, however. 

Fig. [3] also shows the fraction of galaxies at each stellar 
mass which are satellites according to our standard assump- 
tion that M» is a monotonic function of Mhaio (the dashed 
curve which is labelled on the right-hand axis). This frac- 
tion maximises at about 40% for the smallest galaxies con- 
sidered in the plot. Roughly 25% of galaxies of Milky Way 
mass are satellites, but less than 10% of galaxies with M* > 
10 11,5 Mq. (These fractions are consist ent with those in- 
ferred using quite different arguments bv lMandelbaum et alj 
(2006)). These satellite fractions account for at least partly 
the dispersion in the Mhaio - M* relation. 

3.3 Halo masses for galaxies in the Local Group 

Within the Local Group a variety of dynamical tracers are 
available which can provide halo mass estimates for indi- 
vidual galaxies. Clearly, it is of interest to see how these 
compare with the halo masses we infer from our abun- 
dance matching argument. A significant obstacle to carrying 
through this programme is the difficulty of estimating stellar 
masses for comparison to the more distant SDSS galaxies. 
The large angular size of objects like M31, M33 and the 
Magellanic Clouds makes it difficult to obtain accurate inte- 
grated photometry, while our position within the Milky Way 
makes it particularly difficult to infer our own Galaxy's total 
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Name M* M halo 80% Confidence Ref. 

[M Q ] [Mq] Interval [M ] 



M31 


6.98 


X 


10 10 


2.96 


X 


10 12 


[1.03, 7.32] 


X 


10 12 


(1) 


Milky Way 


5.5 x 10 10 


1.99 


X 


10 12 


[0.80,4.74] 


X 


10 12 


(2) 


M33 


2.84 


X 


10 9 


1.74 


X 


10 11 


[1.30,2.23] 


X 


10 11 


(1) 


LMC 


1.30 


X 


10 9 


1.21 


X 


10 11 


[0.93, 1.56] 


X 


10 11 


(1) 


M32 


1.24 


X 


10 9 


1.19 


X 


10 11 


[0.91, 1.53] 


X 


10 11 


(1) 


NGC 205 


9.29 


X 


10 s 


1.05 


X 


10 11 


[0.81, 1.34] 


X 


10 11 


(1) 


SMC 


2.63 


X 


10* 


6.39 


X 


10 10 


[5.23,7.91] 


X 


10 10 


(1) 




5.95 


X 


10 8 


8.69 


X 


10 10 


[0.69, 1.10] 


X 


10 11 


(3) 



Table 1. Stellar and halo masses for selected luminous Local Group galaxies. Stellar masses are computed from B-V colors and V-band 
magnitudes using the fit in table 1 of Bell & de Jong for all galaxies except the Milky Way, for which M, is taken directly from Flynn et 
al. 2006. The colors and magnitudes are taken from the references in column 5. The halo masses in column 3 come from our abundance 
matching results (Fig. 2). while the values in column 4 give the 10% and 90% values of Mhalo assuming cri og M r = 0-2. References: (1) 
de Vaucouleurs et al. 1991 (RC3); (2) Flynn ct al. 2006; (3) Bothun & Thompson 1988. 



stellar mass. In Table 1 we list stellar masses for the brighter 
galaxies of the Local Group, together with the source from 
which we obtained them. Where possible, we have chosen 
estimates based on similar techniques and assumptions to 
those we use for more distant systems. 

For given stellar mass, the M* - Mhalo relation of Fig. [2] 
predicts a unique value of Mhalo- We list this for each galaxy 
in column 3 of Table 1. If this relation does, in fact, have 
some dispersion, then a range of halo masses is consistent 
with any given stellar mass. As noted above, the dispersion 
in log M» at given Mhalo cannot exceed 0.2 if the high-mass 
tail of the SDSS/DR7 mass function is to be reproduced. 
Semi-analytic simulations of galaxy formation within the MS 
(e.g. DLB07) suggest a dispersion which is indeed roughly 
independent of Mhalo, but is somewhat smaller, ~ 0.17, at 
least for systems with Mhalo > 10 Mq. At very low masses, 
the models suggest a rather larger dispe r sion, ~ 0.3 for stel- 
lar masses around 10 s M Q . iMore et alJ (2009) infer a very 
similar dispersion for the luminosity of central galaxies in 
halos of given mass from their SDSS data on satellite statis- 
tics. In column 4 of Table 1 we give the 10% and 90% points 
of the distribution of halo mass predicted, given the stellar 
mass of each galaxy, for a model with the maximal allowed 
dispersion (0.2) in logM, at given Mhalo- 

For the M i lky W ay, the recent stellar mass estimate 
of EH nn et al.l l|2006l ) is quite similar to the old value of 
~ 6 x 10 10 M Q found by 1980's models of Galactic structure. 
For their estimated stellar mass of 5.5 x 10 10 M©, we pre- 
dict the halo mass of the Milky Way to be 2.0 x 10 12 M Q . 
In a model with the maximally allowed dispersion, the up- 
per and lower 10% points of the predicted halo mass dis- 
tribution are at 4.7 X 10 12 Af Q and 0.80 X 1O 12 M . Most 
independent estimates of the halo mass of the Milky Way 
have come from escape velocity arguments, or from Jeans 
modelling of the radial density and velocity dispersion pro- 
files of distant halo tracer populations, e.g. red giants, 
blue horizontal branch stars, globular clusters, or satel- 
lite galaxies. Recent studies have typically come up with 
halo masses in the range 1 to 2 x 10 12 Mp,. Estimates, 
in un its of 10 12 M o , inc lude 1.9 + ? j JWilkinson fc Evans] 
1 19991 ). 2.5+° o or 1 - 8 -o.7 (jSakamoto et al.ll2003t ) depending 
on w hether or not Leo I is included (see also. iLi fc W hite 
|200S| , for the former case), 1.42 + J;^ (|Smith fc et al.ll2007l ). 



0.5 - 1. 5 jBattaglia et al.l [20051 ; iDehnen et al]|2006l ). and 
l.Oto.l (|Xue et al.ll2008l ). These are consistent with the val- 
ues implied by our abundance-matching, though some are 
at the lower end of the permitted range even when the M* 
- Mhalo relation is allowed to have maximal scatter. 

M31 appears to have a larger stellar mass than the 
Milky Way, consistent with its larger maximum rotation ve- 
locity, and this translates into a larger inferred halo mass 
Mhalo = 3.0 x 1O 12 M0. It is interesting that the sum of the 
halo masses estimated for M31 and the Galaxy from our 
"zero scatter" abundance matching is close to the best esti- 
mate of the same qua ntity (M ha i ( MW) + M ha lo(M31) = 
5.3 x 1O 12 M ) which ILi fc White! (|2008i ) obtained from a 
ACDM-calibrated Timing Argument applied to the rela- 
tive orbit of the two galaxies. Note, however, that if we 
allow maximal scatter in the M» - Mhalo relation, then 
the Milky Way's halo mass could be as low as the val- 
ues found in other recent MW analyses, and the sum 
Mhaio(MW) + M ha i (M31 ) would still not vi olate the 90% 
confidence range quoted bvlLi fc White! \ 200ct ). In this case, 
the Milky Way's halo would of course, be substantially less 
massive than those of typical galaxies of similar stellar mass. 

The other Local Group Galaxies listed in Table 1 are 
all predicted to have (maximum past) halo masses at least 
a factor of 10 smaller than those of the two giants. As a 
result, they are likely to have caused relatively little pertur- 
bation to the orbital dynamics of the main binary system. 
The brightest of the satellites are nevertheless predicted to 
have sufficiently massive halos that dynamical friction may 
have modified their orbits. In addition, all of the galaxies 
show evidence for tidal truncation (M32), tidal distortion 
(NGC205, M33), or associated tidal streams (LMC, SMC, 
M33), so it is likely that their current halo masses are smaller 
than the maximum values quoted in Table 1. 



3.4 Stellar mass autocorrelations 

In addition to esti mating the stellar m ass function of galax- 
ies for SDSS/DR7. lLi fc White! (dHI) also studied the clus- 
tering of stellar mass using the same galaxy sample. This 
was quantified by the projected autocorrelation function of 
stellar mass, w*(r p ). On scales larger than individual galax- 
ies Wp(r p ) can be estimated with high accuracy over about 
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Figure 4. The projected stellar mass autocorrelation function in 
the SDSS/DR7 as measured by Li & White (2009) is plotted as 
triangles with error bars and is compared to the result obtained 
if z = (sub)halos in the MS (dashed line) and the MS-II (solid 
line) are populated with galaxies according to the M„ - M ha i a 
relation of Fig. 2. 



Figure 5. Galaxy formation efficiency as a function of halo mass. 
The black curve indicates the values required if a ACDM universe 
is to fit the observed SDSS/DR7 stellar mass function. Coloured 
symbols show the values found for a large number of recent simu- 
lations of the formation of individual galaxies from ACDM initial 
conditions. Different colours correspond to simulations by differ- 
ent authors as noted. The simulation results vary widely, but the 
great majority lock too many baryons into stars to be viable mod- 
els for the bulk of the real galaxy population. 



three orders of magnitude in r v and is rema r kably well de- 
scribed by a single power law. iLi fc White! (I2009T ) showed 
that this behaviour is approximately, but not perfectly re- 
produced by existing galaxy formation simulations. 

In Figure [4] we show the predictions for Wp(r p ) which 
result if (sub)halos in the MS and MS-II (dashed and solid 
lines, respectively) are populated with galaxies according to 
the M, - Mhaio relation of Fig. [2] We compa r e the se esti- 
mates to the SDSS/DR7 results of lLi fc Whitel (120091 ). Over 
the range 20 kpc< r v <20 Mpc, the MS-II prediction is in 
good agreement with the SDSS data. The MS prediction 
converges to the MS-II on scales larger than ~2 Mpc, but 
is significantly too low on smaller scales, becoming roughly 
constant for r p < 100 kpc. This reflects the lower resolu- 
tion of the MS. As noted above, it underpredicts (sub)halo 
abundances for Mhaio < lO 12 A/0 because many of these cor- 
respond to satellite subhalos which have been stripped to 
masses below the MS resolution limit. The objects missed 
are primarily in the inner regions of massive halos, so their 
absence results in a depression of small-scale clustering. In 
semi-analytic galaxy formation simulations based on the 
MS, this effect is addressed by explicitly following "orphan" 
galaxies from the time their subhalos disrupt until the time 
that the code determines that they should t hemselves dis- 
rupt or merge into the cen tral galaxy (see, e.g. lSpringel et ail 
120051 ; ICroton et al.|[20ol ). This effect is negligible in MS- 
II since, as noted above, the subhalo samples are essen- 
tially complete down to Mhaio values that are small enough 
(~ 1O 1O M0) that their galaxies account for almost all stars. 
Thus the excellent agreement between MS-II and the SDSS 
data provides a powerful consistency check on the general 
framework explored in this paper. 



3.5 Galaxy Formation Efficiency 

Given the relation between halo mass and stellar mass, we 
can define a galaxy formation efficiency as the fraction of 
all baryons nominally associated with the halo (calculated 
as the universal baryon fraction times the halo mass) which 
are locked into stars. Thus, 

Efficiency = — x — = 0.17 x — . (4) 

Mhalo lib Mhalo 

We show this galaxy formation efficiency as a function 
of dark matter halo mass in Fig. [5] It peaks at around 
20% in halos with Mhaio ~ 6 x 1O 11 M0, somewhat less 
than the halo mass of the Milky Way. Similar numbers 
have previously been deri ved f rom analogous argum ents by 
iMandelbaum et al l (|2006h and lBaldrv et alj (|2008l ). among 
others. These low efficiencies must be matched by galaxy 
formation simulations if these are to provide a realistic de- 
scription of the formation of real galaxies. In fact, how- 
ever, as shown by the coloured symbols in Fig. [5] most 
recent simulations of the formation of galaxies of Milky 
Way mass convert 25% - 60% of the available baryons 
into stars (Scannapieco et al. 1 120091 : iGovernato et alj l2007l ; 
lOkamoto et al.ll2005l ). The efficiency in I Abadi et all (|2003l ) 
is even higher, due to these authors' neglect of SN feed- 
back. Cyan dots sho w results of a survey of bary onic physics 
parameter space bv lPiontek fc Steinmetj (|2009l ). Several of 
their models do show formation efficiencies as low as re- 
quired to match the SDSS stellar mass function in a ACDM 
universe, but the typical value is around 35%, almost twice 
as large as required. 
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Figure 6. The stellar mass "Tully-Fisher" relation. Small black 
dots plot the stellar masses predicted for individual MS-II halos 
by the M* — M^aio relation of Fig. 2 against their current maxi- 
mum circular velocity. Big symbols are based on estimates of the 
stellar mass of real galaxies from I band (crosses), K band (filled 
circles), and B and R band (open circles) photo metry and on di- 
rect m easures of their maximum rotation velocity iBell &: de Jong 
l200llV 



Galaxy formation efficiency drops very rapidly towards 
both higher and lower mass. In galaxy groups of mass 
1O 13 M0, only 6% of the total baryons can condense to the 
centre and form stars. This reduction in efficie ncy may per- 
haps r eflect the effects of feedback from AGN ()Croton et al.l 
l2006t iBower et all [2006). In halos of mass around 4 x 
1O 1O M0, around 1% of the available baryons have been con- 
verted into_stars. Here, following the original suggestion of 
iLarsonl (|l974l ^ . SN feedback is believed to be responsible 
for the low efficiency, since it can expel gas effectively from 
such shallow potential wells. In the smallest systems, reion- 
ization may also pla y a role in supp r essing condensation 
and star formation (jEfstathioul [l992l ; H cnson et all 120021; 
ISawala et ai1l2009l ). 



3.6 The stellar mass "Tully-Fisher" relation 

A long standing problem in ACDM cosmology has 
been to reproduce simultaneously the galaxy lumi- 
nosity function and the zero-point of the Tully-Fisher 



relation |Kauffmann et al. 


19931 ICole et al.l Il994l: 


iNavarro & Steinmetzl 200d: 


Cole et all 200Ch. The abun- 


dance matching method usee 


; here reproduces the observed 



stellar mass function automatically. To establish a link 
between stellar mass and circular velocity, we use the M* 
- Mhaio relation to assign a stellar mass to the central 
galaxy of each dark matter (sub)halo in the MS-II. The 
maximum circular velocity of these subhalos is tabulated 
in the simulation database, and we take this as a proxy for 
the maximum rotation velocity of the galaxy. The resulting 
"Tully-Fisher" relation is shown in Fig. El Black dots are 
our predictions for central galaxies. There is a tight relation 



between stellar mass and halo maximum circular velocity 
which can be described approximately by a double power 
law. The bend corresponds to the turn-over point in the 
M* - Mhaio relation (Fig. [2}, and to the point where 
galaxy formation efficiency reaches its maximum. Large 
symbols show data fo r observed spiral galaxies taken from 
iBell fc de Jond (|200lh who adopted a scaled-down Salpeter 
Initial Mass Function (IMF) when deriving their stellar 
masses. This gives values which are higher by 0.15 dex than 
those used here, which assume a Chabrier IMF. We shift 
the observational data downwards by 0.15 dex in order to 
compensate for this difference. 

The ACDM model predicts circular velocities which are 
similar to or lower than those observed over the full stel- 
lar mass range. Moderate differences in this sense are ex- 
pected, since the simulations do not a ccount for the gravity 
of the baryons (see the discussion in INavarro fc Steinmetzl 
(2000)). In the region 2.0 < log V c < 2.2 where spiral galax- 
ies dominate the mass functions, the predicted circular ve- 
locity at each stellar mass is lower than in the observations 
by about 25%. This is plausible, giv en results from detailed 
simulations of spi ral formation (e.g. iGustafsson et al" 1 120061 : 
lAbadi et aLll2009T ). These simulations show that galaxy con- 
densation leads to a compression of the inner dark halo 
which is similar to but somewhat weaker th an that predicted 
by simple adiabatic contraction models ([Barnes fc White] 



1 19841 : iBlumenthal et al.ll 19861 : iGnedin et alJl2004h . The com- 
bined effect of the baryonic galaxy and the compressed 
dark halo is an increase in maximum circular velocity which 
may be of the order we require. Such an enhancement is 
already included ap proximately in many disk galax y for- 
mation mode ls (e.g. iMo et al.lll998l: ICole et al.| 200ft) . Re- 
cent studies (jGnedin et al.l 120071 : iDutton et all 120071 ') sug- 
gest that still larger circular velocity enhancements may 
be produced, in which case the galaxy luminosity function 
and Tully-Fisher relation cannot be matched simultaneously, 
but this depends on the details of galaxy formation and 
assembly and goes beyond the issues we can discuss here. 
At higher masses log V c > 2.2, the difference between pre- 
diction and observation is smaller, though the model still 
predicts a slightly lower circular velocity for a given stel- 
lar mass. This again corresponds well to observation since 
strong gravitational lensing studies suggest that the circu- 
lar velocity is constant at close to the halo maximum value 
throughout th e central regions of t hese higher mass galaxy- 
halo systems |Gavazzi et alj 120071 ). At the low-mass end, 
log V c < 2.0, dark matter dominates the gravity throughout 
observed galaxies and the model prediction matches observa- 
tion rather well. In summary, at given stellar mass, circular 
velocities are expected to be higher than represented by the 
small dots in Fig. HJ] because of the gravitational effects of 
the baryons. The observed Tully-Fisher relation could nev- 
ertheless be reproduced if the increase in circular velocity 
is typically ~20%, a value which is not implausible given 
current simulation data. 



4 CONCLUSIONS 

We have related the stellar mass of a galaxy to the dark 
matter mass of its halo by adopting a ACDM cosmology and 
assuming that the stellar mass of a galaxy is a scatter-free 
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and monotonic function of the maximum mass ever attained 
by its halo. By combining the MS and the MS-II we are 
able to derive accurate halo/subhalo abundances for maxi- 
mum masses spanning the entire relevant range, 1O 1O M up 
to 1O 15 M . The SDSS/DR7 provides precise galaxy abun- 
dances spanning almost 4 orders of magnitude in stellar 
mass. By comparing the two we have linked stellar mass 
to halo mass with high formal accuracy over this full stellar 
mass range without any further assumptions about galaxy 
formation physics or halo evolution. 

The ratio of stellar mass to halo mass maximises at 
about 3% for galaxies somewhat fainter than L* , and de- 
creases rapidly towards both higher and lower mass. Except 
possibly for the most massive galaxies, halo masses derived 
from this abundance matching argument agree well with 
those estimated from weak gravitational lensing of back- 
ground galaxies and from a stacking analysis of SDSS satel- 
lite dynamics. We investigated whether the discrepancy at 
high mass might reflect a dispersion in the stellar mass of 
central galaxies at fixed halo mass, and found that the max- 
imal allowed dispersion of 0.2 in log M* at given M^aio leads 
to somewhat better agreement between our prediction and 
observations. 

Using our M* - Mhaio relation we have predicted halo 
masses for a number of the more massive Local Group galax- 
ies. For the Milky Way the inferred halo mass is around 
2 x 10 12 Mq, consistent with most recent estimates from the 
dynamics of halo tracer populations. The inferred halo mass 
for M31 is larger, around 3 x 1O 12 M0. The sum of the two is 
in excellent agreement with the value found when the Tim- 
ing Argument is applied to the relative motion of the two 
giants. The halo masses of the brighter Local Group satel- 
lites are all inferred to be less than 10% of those of the giants, 
and so should have only rather minor effects on Local Group 
dynamics. 

Galaxy formation efficiency peaks at ~ 20% in halos 
slightly less massive than the hosts of ~ L* galaxies. It 
drops rapidly at both higher and lower mass. Similar values 
have been derived previously by others from weak lensing 
and abundance matching studies. Comparison with recent 
hydrodynamic simulations of galaxy formation shows that 
most simulations have conversion efficiencies which are too 
high for them to be viable models for the bulk of the real 
galaxy population. 

The stellar mass - halo mass relation which we derived 
from our abundance matching argument can be used to pop- 
ulate halos in the Millennium and Millennium-II simula- 
tions. For the latter, the implied spatial clustering of stellar 
mass is in remarkably good agreement with a direct and 
precise measurement based on the full SDSS/DR7 dataset. 
Clustering is underpredicted on scales below a Mpc in the 
Millennium Simulation because subhalos corresponding to 
relatively massive satellite galaxies on tightly bound orbits 
are often missed at MS resolution. This effect is mitigated in 
galaxy formation simulations based on the MS by following 
"orphan" galaxies after their subhalos have disrupted. 

The maximum circular velocity of each subhalo in the 
MS-II can be identified with the maximum rotation veloc- 
ity of the central galaxy assigned to it by the above pro- 
cedure. This results in a stellar mass "Tully-Fisher" re- 
lation which we studied over the rotation velocity range 
1.8 <log(V r c [km/s]) < 2.5. For galaxies like the Milky Way, 



this model predicts circular velocities which are about 20% 
lower than observed. This is roughly consistent with the dif- 
ference expected due to the neglect of the gravitational ef- 
fects of the stars. At lower masses dark matter dominates 
throughout the galaxies, and our results match the observa- 
tions quite well. Thus, the ACDM cosmology does seem able 
to reproduce observed luminosity functions and Tully-Fisher 
relations simultaneously. 

Although the abundance matching scheme is a powerful 
way to relate galaxies to their dark matter halos, in reality, 
there must be some scatter in the relation which is likely 
to depend on other physical properties of the systems. An 
important source of scatter is the evolution of the stellar 
mass - halo mass relation. The gas fraction of galaxies is 
very likely greater at high redshift than in the local uni- 
verse. More galactic baryons may be in the form of gas than 
in stars during the early stages of galaxy formation. Fig. [3] 
shows that the satellite fraction is higher for smaller stel- 
lar mass, suggesting there may be more scatter in low-mass 
galaxies than at high mass. Additional scatt er may come 
from variations in halo a ssembly history (e.g. ICroton et al.l 
l2007l ; lGao fc White! 12007m . For example, the star formation 
efficiency is higher in merger-induced bursts than in quies- 
cent phases. We can study such effects by comparing ob- 
servations to models for galaxy correlations as a function 
of redshift and of stellar mass and age. Direct and precise 
measurement of the stellar mass function at high redshift 
will also help us study the scatter more quantitatively, be- 
cause the way in which galaxies populate halos at different 
redshifts is tightly coupled to the evolution of the dark halo 
distribution, and so to the merger trees we construct for our 
simulations. 



ACKNOWLEDGEMENTS 

We thank Maiangela Bernardi and Francesco Shankar 
for useful discussion of stellar mass estimates. The MS 
and MS-II simulations used in this paper were carried 
out at the Computing Centre of the Max Planck Soci- 
ety in Garching. The halo data are publicly available at 
|http://www.mpa-garching.mpg.de/millennium/| 



APPENDIX A 

The stellar mass estimates in lLi fc White! (| 2009) were based 
entirely on the SDSS Petrosian magnitudes. For each galaxy 
these are measured for all five bands within the same circular 
aperture, defined to have radius twice the galaxy's measured 
Petrosian radius. This has the advantages that the colours 
then refer to a specific and well defined stellar population, 
and that no extrapolation from the direct measurements is 
involved. On the other hand, some light falls outside this 
aperture so that the Petrosian magnitudes underestimate 
the total luminosities of galaxies by amounts which depend 
on the shapes of their individual surface brightness profiles. 
For example, the SDSS Petrosian magnitudes recover almost 
all the flux from objects with exponential profiles, but only 
about 80% of the flux from objects with de Vaucou leurs 
profiles. In the procedure of iBlanton fc Roweisl (|2007T l the 
Petrosian colours of each galaxy are fit to stellar population 
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synthesis templates shifted to its known r edshift. Fo r an a s- 
sumed initial mass function (here that of IChabrierl ([2003)) 
this gives a stellar mass-to-light ratio for the galaxy. Mul- 
tiplying by the Petrosian luminosity then gives the stellar 
mass estimate. 

In this paper we are concerned to describe the rela- 
tion between the total stellar mass of galaxies and their 
halo mass. In addition to Petrosian magnitudes, the SDSS 
databases give "model" magnitudes for each galaxy. Expo- 
nential and de Vaucouleurs models are convolved with the 
point-spread-function and fit to the r-band aperture pho- 
tometry. The total magnitude corresponding to the best fit 
of the two is then defined as the r-band model magnitude of 
the object. This clearly involves extrapolation from the mea- 
surements, and in individual cases the errors can be signif- 
icant. Nevertheless, numerous tests suggest that in general 
this produces unbiased and reasonably robust estimates of 
total luminositjQ. We use these model magn i tudes to correct 
the "Petrosian masses" used bv lLi fc White! (|2009h to "total 
masses" by multiplying the former by the ration of the fluxes 
corresponding to the r-band model and Petrosian magni- 
tudes. Note that this procedure ensures that the colours fit 
to the population sysntehsis templates do indeed refer to a 
well defined stellar population which encompasses the bulk 
of each galaxy's stars. This would not be the case if we had 
instead used "model" colours directly. 

In fig. Al we compare the "total" stellar mass func- 
tion (red) obtained in t his way to the "Pe trosian" stellar 
mass function (black) of iLi fc White] (|2009l ). We only plot 
masses above 1O 1O M since the two functions are indistin- 
guishable at lower mass. It is interesting that the differences 
are largest at the highest masses, reflecting the fact that 
massive galaxies tend to have de Vaucouleurs rather than 
exponential profiles. Indeed, the highest mass systems are 
often cD's with extended envelopes which rise above a de 
Vaucouleurs fit to their high surface brightness regions. The 
stellar masses of such systems will still be systematically 
underestimated by our "total" masses. 

As shown in Fig.Al, except at the highest masses, the 
difference between the two mass functions is quite well rep- 
resented by a shift in mass of 0.04 dex. A slightly better rep- 
resentation of the "total" stellar mass function is obtained 
by modifying the model given in ILi fc White! (|2009h to have 
the parameters listed in Table 2. 
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Table 2. Parameters of a triple Schechter function fit to the "total" stellar mass function. 
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